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Abstract 


In the current study the flood risk in the Ojan Chay basin, which is situated in the province 
of East Azarbaijan, Bostan-Abad, was assessed. 16 spatial variables influencing the 
probability of flooding were taken into account, integrated, and examined using fuzzy logic 
within the geographic information system framework in order to extract the risk map. The 
application of the fuzzy gamma operator to the overlapping thematic maps revealed that 
about 4.8% of the basin's area falls into the high risk category. A considerable portion of 
these areas match the basin's major river beds and the surrounding terrain. The class with 
high risk is also located mostly in the vicinity of the mentioned class and includes more 
than 16.6% of the area of the study area. Floods with a ten-year or longer return period may 
damage these areas. In the high risk class, there are roughly 2.3% of residential areas with 
an area of 11.1 hectares. However approximately 42.9% with an area of 35.9 hectares are 
located in the high risk class. The HEC-HMS model was employed in the second phase of 
the study to simulate rainfall-runoff and pinpoint sub-basins with high potential for runoff 
generation. The findings demonstrated that geomorphometric features and land cover have 
a significant impact on the flood hydrograph of the sub-basins and the basin's output. High 
slope, low permeability, a high concentration of rocky outcrops, impermeable surfaces, and 
the least amount of protective cover are the characteristics of sub-basins that exhibit high 
peaks. The upstream sub-basins with the highest peak discharge are 1, 3, 7, 11, and 12. 
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حکىده 


در پژوهش حاضر خطر سیلاب در سطح حوضه اوجان‌چای» واقع در شهرستان بستان‌آباد استان 
آذربایجان شرقی ارزیابی شد. برای تهیه نقشه خطر تعداد VV‏ متغیر مکانی تأثیرگذار بر وقوع سیلاب 
موردتوجه قرار گرفته و در بستر سیستم اطلاعات جغرافیایی با استفاده از منطق فازی ترکیب و تحلیل 
شدند. نتایج همپوشانی با کاربست اپراتور GUIS‏ فازی OUS‏ داد بالغ بر 4/۸ درصد از مساحت حوضه 
در کلاس با خطر بسیار زیاد قرار گرفته است. بخش قابل‌توجهی از این پهنه‌ها منطبق بر بستر رودخانه- 
cle‏ اصلی حوضه و اراضی پیرامون آن‌هاست. کلاس با حطرپذیری زیاد نیز عمدتاً در مجاورت کلاس 
مذکور قرار گرفته و بالغ بر ۱۹/۹ درصد مساحت حوضه مطالعاتی را شامل می‌شود. این پهنه‌ها می‌توانند 
توسط سیلاب‌های با دوره بازگشت ۱۰ سال و بالاتر تحت تأثیر قرار گیرند. حدود ۳/۲ درصد از مناطق 
مسکونی با مساحتی بالغ ۱۱/۱ هکتار در کلاس با خطرپذیری بسیار زیاد و حدود ۶۲/۹ درصد با مساحتی 
بالغ بر ۳۵/۹ هکتار در کلاس با حطرپذیری زیاد واقع شده‌اند. در بخش دوم تحقیق jl‏ مدل HEC-HMS‏ 
به‌منظور شبیه‌سازی بارش - رواناب و شناسایی زیرحوضه‌های با پتانسیل تولید رواناب بالا استفاده شد. 
نتایج نشان داد هیدرو گراف سیلاب زیرحوضه‌ها و خروجی حوضه تا حد زیادی تحت تأثیر ویژگی‌های 
ژئومورفومتری و پوشش زمین است. پیک‌های بالا در زیرحوضه‌های با کمترین میزان پوشش حفاظتی. 
شیب obo‏ نفوذپذیری کم فراوانی برون‌زدهای سنگی و سطوح غیرقابل نفوذ مشاهده می‌شود. می‌توان 
به زیرحوضه‌های بالادست ازجمله ۱. ۳ ۷ ۱ و ۱۲ اشاره کرد که از بالاترین ميزان دبی پیک 
برخوردارند. 

کلیدواژه‌ها: ارزیابی خطر سیلاب نقشه‌سازی خطر سیلاب. سیستم اطلاعات جغرافیایی (GIS)‏ 
« مدل ye HEC-HMS‏ 4 اوجان‌چای» شمال‌غرب ایران 
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مقدمه 


طی ۲۰ سال گذشته. فراوانی مخاطراتی مانند زمین‌لرزه‌هاء» سونامی‌هاء فوران‌های آتشفشانی و حرکات توده‌ای 
تقریباً ثابت باقی مانده است» ولی مخاطرات مربوط به اقلیم (عمدتاً سیلاب‌ها و طوفان‌ها) به‌طور قابل‌توجهی 
افزایش lA‏ کرده است (Centre for Research on the Epidemiology of Disasters, UN International‏ 
(Strategy for Disaster Reduction (CRED, UNISDR) (2015)‏ روند افزایشی رویدادهای اقلیمی حدی 
که ناشی از گرمایش جهانی می‌باشد یکی از مهم‌ترین مسائل جهانی به‌شمار می‌رود که نوع بشر با آن روبروست 
(Chen, et al., 2015; Sui, Lang, & Jiang, 2018)‏ در «e‏ رخدادهای اقلیمی حدی ازجمله سیل و 
خشکسالی بیش از ۷۰ درصد از کل خسارت‌های اقتصادی ناشی از بلایای طبیعی را به خود احتصاص داده‌اند 
(Toya & Skidmore, 2007)‏ در این رابطه. سیل یکی از مخرب‌ترین بلایای طبیعی به‌شمار می‌رود که باعث 
خسارات جانی و مالی قابل‌توجهی در سراسر جهان می‌شود. تا جایی که در سال‌های اخیر حدود ۶۰ درصد از 
کل خسارات اقتصادی ناشی از بلایای طبیعی مربوط به وقوع سیلاب بوده است (یاریان! و همکاران ۲۰۲۰ 
بدین ترتیب. سیل در طی صدها سال گذشته اثرات منفی اقتصادی- اجتماعی در اغلب کشورها برجای گذاشته 
است. در مواجهه با چنین تهدیدی» مطالعه علمی پارامترهای مثر بر سیل و یافتن راه‌حل‌هایی برای پیش‌بینی 
و به حداقل رساندن خسارات ol‏ بر Ole‏ و مال انسان‌ها از اهمیت بالایی برخوردار است. بدین ترتیب. برای 
پیشگیری و کاهش خسارات ناشی از سیل» ارزیابی و تهیه نقشه پهنه‌بندی خطر سیلاب ضروری است (Dung,‏ 
Long, Tran, & Tuyet, 2021)‏ 


TEC‏ کامل de^‏ ممکن نیست؛ ولی ارزیابی و Ce pte‏ سیلاب‌های محتمل آینده امکان‌پذیر است. در سه 
متعاقب آن» ساخت جاده‌ها یا سایر سازه‌ها در مسیر آب. فرونشست به دلیل پمپاژ بیش‌ازحد و درنهایت تغییر 
یا انحراف مسیر Ob m‏ آب )2011 (Kourgialas & Karatzas,‏ روند فعلی و سناریوهای آتی سیل به اطلاعات 
مکانی و زمانی دقیق Gls‏ دارد. از سوی دیگر» شناسایی مناطق آسیب‌پذیر نسبت به سیل و اتخاذ اقدامات 
مقابله‌ای مناسب می تواند خسارات ناشی از of‏ را به‌طور قابل‌ملاحظه‌ای کاهش دهد. علاوه بر این» پهنه‌بندی 
اقدامات کاهش خطر سیل دارد )2021 Kc lea (Allafta & Opp,‏ پیشگیری از سیل فرایندی پیچیده 


است. اما وقوع OT‏ را می‌توان با استفاده از روش‌ها و اقدامات مناسب پیش‌بینی و تا حدی کنترل کرد. در این 
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رابطه» دسترسی به اطلاعات دقیق و به‌هنگام برای جلوگیری از سیل یا حداقل کاهش اثرات مخرب آن نقش 
کلیدی دارد. یکی از این اطلاعات کلیدی, نقشه پتانسیل سیل است (یاریان و همکاران ۲۰۲۰). با توجه به 
موارد Se‏ در پژوهش حاضر خطر وقوع سیلاب در سطح حوضه آبریز اوجان‌چای در شهرستان بستان‌آباد 
استان آذربایجان شرقی مورد تحلیل فضایی قرار گرفت. حوضه مطالعاتی یکی از حوضه‌های سیل‌خیز منطقه 
می‌باشد که در این زمینه می‌توان به سیلاب فروردین‌ماه سال ۱۳۹۲ اشاره کرد. اختلاف ارتفاع obj‏ درجه 
ناهمواری ML‏ شیب زیاد حوضه وجود سازندهای زمین‌شناسی نفوذناپذیر یا با نفوذپذیری کم اثرگذاری 
توده‌هواهای متعدد و دریافت بارش‌های همرفتی- رگباری. مشارکت ذوب برف در دبی رودخانه (مخصوصاً 
در اواخر فصل زمستان و اوایل فصل Ole‏ تغییر کاربری اراضی و تخریب پوشش گیاهی طبیعی حوضه وجود 
سکونتگاه‌های متعدد در مجاورت آبراهه‌های اصلی. دست‌اندازی و تجاوز به حریم آبراهه‌های اصلی حوضه و 
عوامل متعدد دیگر باعث افزایش مخاطره سیلاب در سطح حوضه آبریز اوجان‌چای شده است. بدین ترتیب 
ضرورت دارد که خطر سیلاب در سطح این حوضه مورد ارزیابی و تجزیه‌وتحلیل قرار گیرد. 

در ادامه برخحی از پژوهش‌های مرتبط با ارزیابی خطر سیلاب مورد اشاره قرار گرفته است. (Ghanavati,‏ 
Karam, & Aghaalikhani, (2013)‏ در پژوهشی سعی کردند تا پهنه‌های با بیشترین خطر وقوع سیلاب در 
حوضه فرحزاد را با استفاده از مدل فازی تعیین کنند. در این رابطه. لایه‌های مربوط به شیب انحنا پروفیل و 
پلانیمتریک» ارتفاع, تراکم زهکشی. فاصله از آبراهه» لندفرم‌ها و کاربری اراضی مورد استفاده قرار گرفت. در 
نقشه نهایی» پهنه‌های با Ja‏ بسیار زیاد در پایین‌دست حوضه و منطبق بر دره hel‏ فرحزاد هستند. 
Moradnezhadi, Jourgholami, & Malekian, (2015)‏ به بررسی کارایی مدل هیدرولوژیک HEC-HMS‏ در 
شبیه‌سازی هیدروگراف حوضه‌های ppl‏ جنگلی خیرود پرداختند. کارایی مدل HEC-HMS‏ پس از واسنجی 
آن برای شبیه‌سازی دبی اوج سیلاب در حوضة مطالعه‌شده مورد تأیید قرار گرفت. ,۲۵0006 Sistani‏ 
Negaresh, & Fotoohi, (2017)‏ مناطق سیل گیر را در حوضه آبریز گابریک در جنوب شرق ایران شناسایی و 
پهنه‌بندی کردند. نتایج حاصل از پهنه‌بندی با استفاده از روش تحلیل سلسله مراتبی OUS‏ داد که ۲۳ درصد از 
مساحت حوضه در گروه مناطق با das‏ بالا و بسیار بالا قرار گرفته است. عواملی همچون شیب و فرسایش- 
پذیری 255 و نفوذپذیری e$‏ از دلایل اين وضعیت می‌باشند. )2022( Ganji, Gharechelou, & Ahmadi,‏ 
خطر سیلاب را در رودخانه گرگانرود در شهرستان آققلا پهنه‌بندی کردند. OLE!‏ ۱۸ معیار موثر بر سیل‌خیزی 
را با استفاده از تحلیل سلسله مراتبی وزن‌دهی کردند و نتایج نشان داد معیارهای توپوگرافی و مورفولوژی در 
ایجاد سیل در حوضه موردبررسی بالاترین تأثیر را دارند. نتایج پهنه‌بندی نیز نشان داد حدود VE‏ درصد از 
مساحت شهرستان آق‌قلا در کلاس he‏ کم و YI‏ درصد آن در پهنه dam‏ خیلی 365 قرار دارند. Mesri‏ 
Alamdari, (2021)‏ خطر سیلاب را در حوضه آبریز قلعه‌چای با استفاده از GIS‏ و مدل HEC-HMS‏ تحلیل 
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کردند. در این تحقیق تعداد هشت متغیر جهت پهنه‌بندی dam‏ سیلاب با استفاده از منطق فازی ترکیب شدند. 
مدل هیدرولوژیکی HEC-HMS‏ نیز جهت شناسایی پهنه‌های با تولید رواناب بالاتر به کار برده شد. نتایج 
حاصل از پهنه‌بندی خطر سیلاب نشان داد که ۷/۷ درصد از مساحت حوضه قلعه‌چای از خطرپذیری بالایی در 
مقابل سیلاب برخوردار است. شبیه‌سازی بارش- رواناب OLS‏ می‌دهد که بارش‌های با دوره بازگشت بالاتر 
می‌تواند دبی‌های سیلابی ایجاد کند. Zare, Azari, & Rezaei Arefi,(2024)‏ مناطق آسیب ‌پذیر نسبت به سیلاب 
را در حوضه آبخیز کشف رود با استفاده از ۵ روش تصمیم‌گیری چندمعیاره شامل تحلیل شبکه» ویکور 
تاپسیس تحلیل سلسله مراتبی و الکتره و با در نظر گرفتن VA‏ شاخص مورفومتری مستخرج از مدل رقومی 
ارتفاع منطقه شناسایی کردند. نتایج پژوهش OUS‏ داد که در بین روش‌های موردبررسیء تحلیل شبکه برای تهیه 
نقشه سیل‌خیزی دارای دقت بالاتری می‌باشد. )2011( «Ali Mohammad, Aslam, & Khan,‏ برای بررسی 
میزان تأثیر تغییر کاربری اراضی بر رواناب سطحی در حوضه لینولاه در پاکستان از مدل بارش رواناب HEC-‏ 
5 استفاده کرد. نتایج نشان داد که در آینده با تغییر کاربری اراضی. رواناب کل بین ۵۱/0 تا ۱۰۰ درصد 
افزايش خواهد پافت. )2011( Nandalal & Ratnayake,‏ در تحقیقی به ارزیابی کارایی رویکرد فازی به‌منظور 
بررسی ریسک سیلاب در سطح حوضه آبریز رودخانه کالوگانگا در سریلانکا پرداختند. گستره سیلاب برای 
بارندگی با دوره باز گشت ۱۰۰ ساله با استفاده از سیستم مدل‌سازی هیدرولوژیک (HEC-HMS)‏ و مدل‌های 
مبتنی بر سیستم تجزیه‌وتحلیل رودخانه HEC-River‏ تعیین شد. مقایسه نتایج به دست آمده نشان داد که منطق 
فازی عدم اطمینان در تعیین سطح خطر, آسیب‌پذیری و سطح خطر را ب‌صورت دقیق‌تری در نظر می‌گیرد. 
Ajin, Krishnamurthy, Jayaprakash, & Vinod, (2013)‏ مخاطره سیل را در حوضه رودخانه واماناپورام در 
هند بر اساس ارزیابی چند معیاره و کاربست تکنیک‌های دورسنجی و سیستم اطلاعات جغرافیایی (GIS)‏ مورد 
ارزیابی قرار دادند و منطقه موردمطالعه را از نظر پتانسیل خطر سیلاب به پنج پهنه مجزا تقسیم کردند. 
Rashetnia, (2016)‏ به ارزیابی آسیب‌پذیری نسبت به سیلاب در مورلندسیتی واقع در منطقه ملبورن استرالیا با 
کاربست Ga»‏ فازی پرداخت. یافته‌های این تحقیق نشان‌دهنده کارایی بالای منطق فازی و سیستم اطلاعات 
جغرافیایی (GIS)‏ برای ارزیابی آسیب‌پذیری نسبت به سیلاب است. این پژوهش نشان داد که مناطق بسیار 
پرجمعیت در طبقات آسیب‌پذیری 265 و بسیار زیاد واقع‌اند. )2020( Tella & Balogun,‏ جهت ارزیابی مکانی 
خطر سیل در ایبادان در نیجریه از ترکیب روش‌های تصمیم‌گیری چندمعیاره و فازی استفاده کردند. برای بهبود 
دقت وزنی معیارهای ایجاد سیل فازی AHP‏ با سیستم اطلاعات جغرافیایی برای محاسبه وزن و فرآیند 
همپوشانی OS‏ شد. نقشه فازی AHP‏ دارای /٩۱‏ سازگاری با مکان‌های تاربخی سیل بود. مروری بر پیشینه 
پژوهش نشان می‌دهد که محققان با استفاده از روش‌های متعددی به ارزیابی خطر سیلاب پرداخته‌اند. این امر 


در طی سال‌های اخیر با تکیه بر تکنولوژی سیستم اطلاعات جغرافیایی صورت گرفته است. در پژوهش حاضر 
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ما dades pet‏ برای ig‏ سا Si sie‏ ی ساب از cont dta‏ اماد 
شده است. ارزیابی سیل‌خیزی يا پتانسیل تولید رواناب نیز متکی بر مدل HEC-HMS‏ در ترکیب با GIS‏ 


illia 


مواد و روش‌ها 
viieni dias‏ 


در پژوهش حاضر خطر وقوع سیلاب در سطح حوضه آبریز اوجان‌چای. واقع در شهرستان بستان‌آباد استان 
آذربایجان شرقی. ارزیابی شده است. از نظر موقعیت جغرافیایی حوضه مطالعاتی در EV‏ درجه و ۲۱ دقیقه و 
asl ۰‏ تا ۶7 درجه و EV‏ دقیقه و £0 asl‏ طول شرقی و ۳۷ درجه و ٤٤‏ دقیقه و ۸ ثانیه تا ۳۷ درجه و ۵۱ 
دقیقه و ۵ انیه عرض شمالی واقع شده است (شکل ۱). مساحت حوضه حدود ۷۲ کیلومترمربع و محیط OF‏ 
حدود ۸۶ کیلومتر می‌باشد. حوضه اوجان‌چای بخش‌هایی از دامنه‌های شمال شرقی توده آتشفشانی سهند را 
زهکشی می کند و از پایکوه‌ها تا نزدیک قله سهند امتداد یافته است. در محدوده حوضه مطالعاتی تعداد ۱۵ 


روستا وجود دارد که اکثراً در مجاورت پا بستر دره‌های اصلی حوضه استقرار یافته‌اند. افزون بر این» دره‌های 
اصلی حوضه آبریز اوجان‌چای کانون فعالیت‌های کشاورزی منطقه می‌باشند. 
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شکل ۱- موقعیت حوضه آبریز اوحان‌جای در شمال غرب ol ul‏ 


۲12.1. The location of the Ojan Chay watershed in the northwest of Iran 
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روش انجام پژوهش 


متغیرهای متعددی بر وقوع سیلاب در یک حوضه آبخیز تأثیر می‌گذارند. مهم‌ترین عوامل و متغیرها عبارتند 
از: عوامل فیزیکی (مانند ارتفا شیب. جهت شیب انحنای دامنه» سنگ‌شناسی» موقعیت توپوگرافی و بارند گی)» 
انسانی (مانند کاربری اراضی و فاصله از جاده) ULL)‏ و همکاران. ۰ با توجه به موارد مذ کور در پژوهش 
حاضر به‌منظور ارزیابی سیل گیری (پهنه‌بندی سیلاب) و سیل‌خیزی (تولید رواناب) حوضه آبریز اوجان‌چای از 
داده‌هایی بدین شرح استفاده به عمل آمد: نقشه‌های توپوگرافی مقیاس ۱:۰۰۰۰ (سازمان جغرافیایی نیروهای 
مسلح) نقشه‌های زمین‌شناسی مقیاس ۱:۱۰۰۰۰۰ (سازمان زمین‌شناسی کشور)؛ نقشه خاک (سازمان آب منطقه- 
ای)» تصاویر مدل رقومی ارتفاع (DEM)‏ مربوط به ماهواره استر" با قدرت تفکیک مکانی ۲۷ متر» تصاویر 
ماهواره‌ای d‏ ۲ (مربوط به اوایل تابستان ٤١١‏ ه. 4( 3 گوگل ارث. داده‌های ایستگاه سینویتیک 
بستان‌آباد و ایستگاه باران‌سنجی دیزناب و داده‌های ایستگاه‌های هیدرومتری بستان‌آباد و دیزناب. ارزیابی مکانی 
سیلاب مستلزم مشارکت عوامل متعددی بوده و می‌بایست به‌صورت ترکیبی و یکپارچه مورد تجزیه‌وتحلیل 
قرار گیرند. بدین منظور در پژوهش حاضر مخاطره سیلاب در سطح حوضه آبریز اوجان‌چای با بهره‌گیری از 
نمودن و ترکیب لایه‌های موضوعی موثر بر سیل‌گیری از رویکرد منطق فازی استفاده شد. همچنین به‌منظور 
شبیه‌سازی بارش- رواناب در سطح حوضه و شناسایی پهنه‌های با بیشترین ميزان تولید رواناب از مدل 
هیدرولوژیکی  HEC-HMS‏ اکستنشن آن در محیط نرم‌افزار ArcGIS‏ یعنی HEC-GeoHMS‏ بهره گرفته شد. 


منطق فازی 

تئوری مجموعه‌های فازی و منطق فازی به‌عنوان نظریه‌ای ریاضیاتی برای مدل‌سازی و صورت‌بندی ریاضی 
ابهام و عدم قطعیت موجود در فرایندهای شناختی انسانی» ابزارهای بسیار کارآمد و مفیدی به شمار می‌رود. 
این نظریه قادر است به بسیاری از مفاهیم» متغیرها و سیستم‌های نادقیق (فازی) و مبهم» صورت‌بندی ریاضی 
بخشیده و زمینه را برای استنتاج» کنترل و تصمیم‌گیری در شرایط عدم اطمینان فراهم آورد Jahanbakhsh Asl,‏ 
Asad, & Akbari, (2016)‏ فرایند فازی‌سازی. عدم قطعیت کلاس‌ها را در خحصوص پدیده‌هایی مشخص 
می‌سازد که دارای مرزهای دقیق نیستند. در فرایند فازی‌سازی, ایده‌آلی برای عضویت در مجموعه تعریف 
می شود. به هر مقداری از پدیده که سبت به هسته مجموعه» محوری تر (مرکزی‌تر) باشد ارزش ۱ احتصاص 


! - Aster 
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بين این دو حد ih gt ga e»‏ در منطقه گذار مجموعه قرار می گیرند Sarvar & Kherizadeh Arouq,(2015)‏ 
تابع عضویت در مدل فازی به‌صورت رابطه (۱) تعریف می‌شود: 

رابطه (۱) (X )} foreach x € X...]‏ ریم , A —(X‏ 
ab is‏ عضویتی است که درجه عضویت X‏ را در مجموعه A‏ تعریف می کند. بعد از فازی‌سازی لایه‌های 
موضوعی نسبت به ترکیب يا روی‌هم گذاری فازی آن‌ها اقدام شد. بسته نرم‌افزاری ArcGIS‏ پنج عملگر GIB‏ 
مختلف را جهت ترکیب فازی لایه‌های رستری در اختیار می‌گذارد این عملگرها به شرح زیرند )2007 (Lee,‏ 
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Ucombination = (Fuzzy algebraic sum)^ * (Fuzzy algebraic product) ^ 


مدل هیدرولوژیکی HEC-HMS‏ 


در این تحقیق جهت شبیه‌سازی بارش - رواناب و بررسی پتانسیل تولید رواناب حوضه اوجان‌چای از روش 
مدل‌سازی هیدرولوژیکی HEC-HMS‏ همراه با الحاقی HEC-GeoHMS‏ استفاده شد. این سیستم مدل‌سازی 
هیدرولوژیک برای شبیه‌سازی فرایندهای بارش- رواناب سیستم‌های حوضه آبریز طراحی شده است. طراحی 
سیستم به نحوی بوده تا در یک طیف گسترده‌ای از مناطق جغرافیایی برای حل مسائل متعدد قابل‌اجرا باشد. 
این گونه مسائل شامل تأمین آب و هیدرولوژی سیلاب حوضه‌های بزرگ رودخانه و رواناب حوضه‌های آبریز 
کوچک شهری یا طبیعی می‌شود. هیدروگراف‌های تولید شده توسط این برنامه Sy pod‏ مستقیم یا eli‏ با 
نرم‌افزارهای So‏ برای مطالعات قابلیت دسترسی آب. زهکشی شهری» پیش‌بینی جریان, اثرات آتی شهرنشینی» 
طراحی سرریز مخزن, کاهش خسارات سیلاب. منطقه‌بندی دشت سیلابی و عملیات سیستم‌ها مورد استفاده 
قرار می گیرند )2016 (Hydrologic Engineering Center,‏ در این HEC-GeoHMS akl)‏ به‌عنوان یک ابزار 
هیدرولوژی زمین- مکانی برای مهندسان و هیدرولوژیست‌ها با تجربه محدود GIS‏ توسعه یافته است. این 
al,‏ نمایان‌سازی اطلاعات مکانی» مستند کردن ویژگی‌های حوضه آبریز» تحلیل مکانی» (ex‏ حدود 
زیرحوضه‌ها و آبراهه‌ها؛ ایجاد ورودی‌ها به مدل‌های هیدرولوژیک و تهیه گزارش‌ها برای کاربران را امکان‌پذیر 
می‌سازد )2013 (Hydrologic Engineering Center,‏ برای مدل‌سازی فرایند بارش - رواناب با استفاده از 
نرم‌افزار HEC-HMS‏ لازم است که سه مدل شامل مدل حوضه مدل هواشناسی و مدل مشخصات کنترل 
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خاک آمریکا (SCS)‏ صورت گرفت. مدل مذکور در کتاب‌های هیدرولوژی به‌تفصیل تشریح شده است؛ بنابراین 
با توجه به حجم زیاد مطالب از بیان مدل مذکور در مقاله حاضر اجتناب شد. 


نتایج و بحث 


پهنه‌بندی daz‏ سیلاب حوضه اوجان‌چای با استفاده از منطق فازی 


برای eg‏ نقشه خطر سیلاب در سطح حوضه اوجان‌چای ابتدا توزیع فضایی مجموعه‌ای از عوامل مژثر بر خطر 
سیلات مورد ارزیابی قرار گرفت. سپس عوامل مذکور با استفاده از aly‏ مختلف در بستر GIS‏ فازی شدند. 
Cha alg)‏ فازی ا انتفاده از lay gl gl‏ مختلت فاوی خر کیب و AAS anal SS‏ 


ارزیابی توزیع فضایی عوامل مؤثر بر خطر سیلاب 


در این بخش از پژوهش توزیع فضایی عوامل مؤثر بر خطر سیلاب در سطح حوضه اوجان‌چای مورد ارزیابی 
قرار گرفته است. در شکل (۲) توزیم فضایی عوامل مربوطه ارائه شده است. 

- ارتفاع: بر اساس محاسبات پیکسل‌پایه» میانگین ارتفاع حوضه اوجان‌چای حدود ۲۲۶۹ متر می‌باشد. همچنین 
اختلاف ارتفاعی بالغ بر ۱۷۹۶ متر در سطح حوضه دیده می‌شود. سطوح مرتفع حوضه از نظر تولید رواناب و 
سیل‌خیزی حائز اهمیت زیادی می‌باشند. رواناب تولیدی به سمت اراضی پست و کم‌ارتفاع حوضه هدایت 
می‌شوند و به دلیل وجود سکونتگاه‌هاه تأسیسات و اراضی کشاورزی در پهنه‌های مذکور می‌توانند باعث 
مخاطراتی شوند؛ بنابراین. سطوح کم‌ارتفاع حوضه از نظر خطر سیل گیری دارای اهمیت بوده و در فرایند 
فازی‌سازی ارزش‌های زیادتری را به خود اختصاص می‌دهند. 

- شیب: متوسط شیب در سطح حوضه اوجان‌چای حدود ۱۷ درصد می‌باشد. در رابطه با پهنه‌بندی. اراضی 
کم‌شیب همواره در معرض پخش سیلاب قرار دارند. در این ably‏ طبقه شیب ۰ تا ۵ درصد از اهمیت بیشتری 
برخوردار است. در حوضه اوجان‌چای بالغ بر ۱۲ درصد مساحت حوضه در طبقه شیب مذکور قرار گرفته 
است. سطوح شیب مذکور از یک طرف در معرض سیلاب‌های رودخانه‌ای می‌باشند و از طرف دیگر می‌توانند 
توسط سیلاب‌های ورقه‌ای (شیت‌فلادها) تحت تأثیر قرار گیرند. از نظر تولید رواناب. عکس موارد مذکور 
Gale‏ است. به‌بیان‌دیگر اراضی پرشیب قسمت‌های میانی و بالادست حوضه از بیشترین میزان تولید رواناب 


برخوردارند. 
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- جهت شیب: اثر مستقیم جهت شیب به‌واسطه تأثیر بر ذوب برف و اثر غیرمستقیم آن به‌واسطه تأثیر بر پوشش 
گیاهی است. در دامنه‌های جنوبی و غربی» برف‌های حوضه با سرعت بالاتری ذوب می‌شوند. این امر در اواخر 
زمستان و اوایل بهار که منطقه بارش‌های همرفتی- رگباری فراوانی را تجربه می کند حائز اهمیت زیادی است. 
تراکم پوشش گیاهی اندک و درنتیجه اثرات محدود حفاظتی آن در مقابله با سیلاب را می‌توان به موارد مذکور 
اضافه نمود. 
- تحدب سطح زمین: شاخص مذکور برای اولین بار توسط )2007( Iwahashi & Pike,‏ مطرح شد و هدف 
از ارائه ol‏ طبقه‌بندی عوارض سطح زمین بود. تحدب سطح زمین معیاری مناسب برای ارزیابی پتانسیل پخش 
یا تجمع آب به‌شمار می‌رود. مقادیر SVL‏ شاخص تحدب زمین معرف سطوح با تحدب زیاد و مقادیر پایین‌تر 
نشان‌دهنده سطوح با تحدب Gal‏ (مقعر) می‌باشد. پهنه‌های با تقعر بالا شامل سطوح با مقادیر ۲۵ تا ۳۱ 
می‌باشند که در حدود ۱۲/۲ درصد از سطح حوضه اوجان‌چای را به خود اختصاص داده‌اند. این پهنه‌ها مشتمل 
بر سطوح مقعر حوضه بوده و به‌طور عمده منطبق بر بستر دره‌ها و سیرک‌های یخچالی بالادست حوضه 
می‌باشند. به دلیل ریخت‌شناسی مقعر» بخش قابل‌توجهی از رواناب‌های تولید شده به سطوح مذکور هدایت 
AP dius‏ 

- طول دامنه: عامل طول دامنه یکی از عوامل کلیدی تأثیرگذار بر حجم جریان و سرعت رواناب به‌شمار 
می‌رود. هر چه طول دامنه بیشتر باشد. سرعت = Ob‏ بالاتر و حجم رواناب بیشتر می‌شود. این امر به‌طور 
قابل توجهی شکل‌گیری سیل را تحت SU‏ قرار می‌دهد )2020 (Dung, Tuyet, Ahmad, Long,‏ در تحقیق 
حاضر برای استخراج aY‏ طول دامنه از رابطه زیر استفاده شد )2011 (Luca, Conforti, & Robustelli,‏ 

رابطه LSF-(fa*Cell size/22.3)^ 0.4*(Sino/0.0896)^ 1.3 (V)‏ 
که در آن: LSF‏ عامل طول دامنه و 18 تجمع جریان می‌باشد. در حوضه مطالعاتی polis‏ طول دامنه از حداقل 
صفر تا حداکثر ٩۰‏ متر متغیر است. بین طول دامنه و وقوع سیلاب رابطه‌ای مستقیم وجود دارد. 

- توان رودخانه: polis‏ بالاتر شاخص مذکور بیانگر این است که جریان رودخانه از توانایی انجام کار 
ژئومورفیک بیشتری برخوردار است. این امر از نظر تولید رسوب. فرسایش و افزایش توان تخریبی سیلاب‌ها 
JE‏ اهمیت می‌باشد. توان رودخانه را می‌توان از Gob‏ فرمول زیر محاسبه نمود & (Luca, Conforti,‏ 
:Robustelli, 2011)‏ 

رابطه )^( SPI-In(As/tanc)‏ 
که در آن: SPI‏ شاخص قدرت جریان» As‏ مساحت ویژه جریان و 6 درجه شیب می‌باشد. در سطح حوضه 
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ارتباط با تغییرپذیری شدید توپوگرافی منطقه می‌باشد. در این رابطه» در سطح حوضه آبریز اوجان‌چای مقادیر 
توان رودخانه از حداقل صفر تا حداکثر ۱۵۰۰۰ وات بر مترمربع متیر است. 

- رطوبت توپوگرافیک: این شاخص یکی از کارآمدترین شاخص‌های مورد استفاده برای ارزیابی خطر سیلاب 
به‌شمار می‌رود. شاخص Se TWH‏ بر این فرض Col‏ که توپوگرافی کنترل کننده حرکت آب در چشم‌انداز 
است. در واقع» شاحص مذکور تأثیر توپوگرافی محلی در فرآیندهای هیدرولوژیکی را به‌صورت کمی OUS‏ 
می‌دهد. شاخحص رطوبت توپوگرافیک از طریق رابطه زير محاسبه می‌شود  (Kumar & Gupta, 2016; Qin,‏ 
al., 2011)‏ : 


(As * Pixel area) 
Tan((S * 1)/180)) 


TWI = Ln (4) رابطه‎ 


که در آن: As‏ مساحت زهکشی دامنه فوقانی در واحد عرض خط کانتور (مترمربع/ متر) و S‏ شیب به درجه 
است. مقادیر این شاخص در سطح حوضه مطالعاتی از حداقل ۳ تا حداکثر ۲۳ متغیر است. رابطه‌ای مستقیم 
بین مقادیر این شاخص و خطرپذیری سیلاب در سطح حوضه برقرار است. در پهنه‌هایی که این شاخص مقادیر 
بالایی دارد احتمال وقوع سیلاب بیشتر است. در این رابطه کلاس با TWI polis‏ بیش از ۱۵ حائز اهمیت 
پیشتری می‌باشد. 

- شاخص همگرایی: شاحص همگرایی پارامتری است که ساختمان ناهمواری را به‌صورت مجموعه‌ای از 
سطوح همگرا و واگرا نمایش می‌دهد. مقادیر منفی این شاخص نشان‌دهنده سطوح همگرا بوده که بیشتر در 
معرض سیلاب قرار می‌گیرند. در این رابطه کلاس با مقادیر همگرایی کمتر از ۲۰- معرف سطوح با همگرایی 
بسیار زیاد هستند. سطوح مذکور بالغ بر Y‏ درصد مساحت حوضه را شامل می‌شوند. این سطوح مشتمل بر 
بستر آبراهه‌ها و اراضی مجاور آن‌ها بوده و همه‌ساله توسط دبی‌های معمول و طغیانی (سیلاب‌های با دوره 
بازگشت Y‏ تا ۳ ساله) تحت تأثیر قرار می‌گيرند. کلاس با مقادیر شاحص همگرایی ۲۰- تا -٤‏ نیز دلالت بر 
سطوح با همگرایی زياد هستند که حدود ۱۳ درصد از مساحت حوضه را به خود اختصاص داده‌اند. با توجه 
به مطالعات میدانی و داده‌های هیدرومتری می‌توان گفت که سطوح مذکور در معرض سیلاب‌های با دوره 
باز گشت Y‏ تا ۲۵ ساله می‌باشند. 

- عمق دره: مقادیر بالاتر این شاخص دلالت بر عمق زیاد دره نسبت به خطالرأس‌ها دارد. این امر به دلایلی 
منجر به افزايش خطرپذیری نسبت به سیلاب می‌شود. در این رابطه می‌توان به مشارکت قابل توجه رواناب‌های 
دامنه‌های مشرف به این دره‌ها در دبی آبراهه‌های اصلی اشاره نمود که حطر وقوع سیلاب را به‌صورت محلی 
افزایش می‌دهد. بدین ترتیب طبقه با مقادیر بالاتر از ۲۳۰ متر از بیشترین میزان خطرپذیری نسبت به سیلاب 


برخوردار اسنت: 


۱۶۰۳ تابستان‎ Y شماره‎ AY نشریه جغرافیا و مخاطرات محیطی» جلد‎ v1 


- تراکم زهکشی: هرچه تراکم زهکشی حوضه‌ای بالاتر باشد تخلیه رواناب‌ها راحت‌تر صورت گرفته و زمان 
تمرکز زیرحوضه‌ها کاهش می‌یابد. این امر منجر به افزایش Jam‏ سیلاب می گردد. مقادیر این شاخص در سطح 
حوضه مطالعاتی از ۰ تا ۵ کیلومتر در کیلومترمربع متغیر است. در کل از رس قلل و ستیغ‌های منطقه در جهت 
cl‏ دامنه‌ها و دره‌ها بر میزان تراکم زهکشی افزوده می‌شود. 

- فاصله از رودخانه: نواحی مجاور رودخانه‌ها بیشتر تحت تأثیر سیلاب‌های رودخانه‌ای قرار می‌گیرند و با 
دور QU‏ او نز E ads lada,‏ واا می کل در این وان جارد ۳۱ درد از elus‏ 
کاربری‌های مسکونی (حدود ۳۰/۲ هکتار) در فاصله صفر تا ۲۰۰ متری از رودخانه‌های حوضه استقرار یافته‌اند. 
پهنه‌های مذکور جزو خطرپذیرترین و آسیب‌پذیرترین پهنه‌های حوضه مطالعاتی نسبت به سیلاب می‌باشند. 
- کاربری اراضی: کاربری‌هایی که باعث افزایش نفوذ آب می‌شوند به‌عنوان کاربری‌های تعدیل کننده سیلاب‌ها 
مطرح می‌شوند. در این رابطه می‌توان به مراتع با پوشش گیاهی خوب اشاره نمود که درصد قابل‌توجهی از آب 
ناشی از بارش را جذب می‌کنند. برعکس» کاربری‌هایی مانند کاربری مسکونی با کاهش نفوذ و افزایش آبدوی 
باعث تبدیل بخش قابل‌توجهی از آب بارش به رواناب سطحی می‌شوند و درنتیجه خطر سیلاب را افزایش 
می‌دهند. 

- پوشش گیاهی: در این زمینه از شاخحص NDVI‏ استفاده شد که از رابطه زیر به دست می‌آید (لئو و ماسون 
(Yard‏ 


NIR — Red (V+) abl 
VI .سس‎ 22 
ADNI NIR + Red 


که در آن: NIR‏ و RED‏ به ترتیب بازتاب‌ها در باند مادون‌قرمز نزدیک و باند ja‏ هستند. بین تراکم پوشش 
گیاهی و خطر سیلاب رابطه‌ای معکوس برقرار است. این hel‏ در فرایند فازی‌سازی لایه NDVI‏ به کار بسته 
شد و بدین ترتیب» پیکسل‌های با مقادیر 22۷1 اندک از ارزش بیشتری برخوردار BAS‏ و برعکس. 

- زمین شناسی: سازندهای زمین‌شناسی حوضه اوجان‌چای را با توجه به میزان نفوذپذیری و تأثیر OT‏ بر وقوع 
سیلاب می‌توان به سه گروه طبقه‌بندی نمود: سازندهای با نفوذپذیری زیاد تا بسیار زیاد: این سازندها به ترتیب 
افزایش نفوذپذیری شامل پادگانه‌های آبرفتی قدیمی (QU)‏ پادگانه‌های آبرفتی Loe‏ (012) و آبرفت‌های عهد 
حاضر (Qa)‏ می‌باشند. در طی کمی‌سازی کدهای ۱ تا ۳ به این سازندها اختصاص cdl‏ سازندهای با 
نفوذپذیری متوسط: این سازندها بر اساس کاهش میزان نفوذپذیری شامل لاهار و جریان‌های گلی Plvc2la)‏ با 
کد dE‏ توف برشی و توف‌های لایه‌ای Et)‏ با کد ۵) و پومیس خاکستر آتشفشانی و کنگلومرا Plvclp)‏ با JS‏ 
1 می‌باشند. سازندهای با نفوذپذیری کم تا بسیار کم: این سازندها روی‌هم‌رفته بالغ بر ۵۸/۱ درصد از مساحت 


حوضه را شامل می‌شوند و عبارتند از: واحد ۳۱۷020 (کنگلومرا با کد (V‏ واحد Eit‏ (ایگنیمبریت- توف با کد 


۱ - Liu & Mason 


سال سیزدهم اسمعیل پور» ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از منطق فازی و مدل هیدرولوژیکی VV  ...‏ 


CA‏ واحد TG‏ (گرانیت الکالی با کد 4 واحد EV‏ (تراکی آندزیت- بازالت با کد ۱۰) و واحد ۲128 (داسیت- 
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شکل ۲- نقشه متغیرهای مورد استفاده در پهنه‌بندی daz‏ سیلاب حوضه آبریز اوجان‌چای 


۲۱2۰2۰ Map of variables used in flood risk zoning of Ojan Chay watershed 


۷۸ نشریه جغرافیا و مخاطرات محیطی» جلد AY‏ شماره ۰۲ تابستان ۱۶۰۳ 


تراکم زهکشی N‏ عمق دره (به متر) 
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سازندهای زمین شناسی کاربری اراضی 
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ادامه شکل ۲- نقشه متغیرهای مورد استفاده در پهنه‌بندی خطر سیلاب حوضه آبریز اوجان‌چای 
Map of variables used in flood risk zoning of Ojan Chay watershed‏ ۲۱2۰2۰ 
فازی‌سازی لایه‌های موضوعی FB‏ بر سیلاب 


فازی‌سازی لایه‌های موضوعی با هدف بی‌بعدسازی و همسان‌سازی لایه‌هاء ارزش گذاری مجدد WAY‏ و رفع 


عدم قطعیت‌های کلاسه‌بندی bad‏ می‌باشد. در این رابطه لایه‌های موضوعی مورد استفاده در پهنه‌بندی خطر 


سال سیزدهم اسمعیل پون ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از منطق فازی و مدل هیدرولوژیکی .. ۰ ۷۹ 


سیلاب را می‌توان به دو گروه طبقه‌بندی نمود: الف) لایه‌هایی که در آن‌ها ارزش‌های بالاتر منجر به افزایش 
خطر سیلاب می‌شوند. توان رودخانه. طول دامنه. عمق دره. رطوبت توپوگرافیک و تراکم زهکشی در این گروه 
قرار می گیرند. درنتیجه. به‌منظور فازی‌سازی Gad‏ مذکور می‌توان از توابعی مانند تابع خطی افزایشی b‏ 
بزرگ فازی استفاده نمود. ب) لایه‌هایی که در آن‌ها ارزش‌های پایین‌تر منجر به افزايش خطر سیلاب می‌شوند. 
ارتفا شیب. تحدب سطح زمین. همگرایی. فاصله از آبراهه و پوشش گیاهی در این گروه قرار می‌گيرند. دو 
تابع فازی خطی کاهشی و نزدیک فازی می‌توانند برای فازی‌سازی لایه‌های مذکور مورد استفاده قرار گیرند. 
سه لایه کاربری اراضی. لیتولوژی و جهت شیب لایه‌های کیفی تحقیق می‌باشند. ازاین‌رو در گام نخست لایه‌های 
موضوعی مذکور از طریق کددهی کمی گردیده و سپس فازی شدند. در شکل (۳) نمونه‌ای از نقشه فازی ارائه 


شده است. 


شکل ۳- نمونه‌ای از نقشه فازی تحقیق (لایه ارتفاع) 
Fig.3. An example of a research fuzzy map (elevation layer)‏ 


ترکیب فازی bal‏ و angi‏ نقشه خطر سیلاب 


در پژوهش حاضر لایه‌های موضوعی فازی با کاربست پنج اپراتور «و» فازی. thy‏ فازی. جمع جبری $38« 
ضرب جبری فازی و فازی گاما ترکیب شدند و نقشه خطر سیلاب حاصل گردید. در پهنه‌بندی با اپراتورهای 
«و» فازی و ضرب فازی اکثریت منطقه در کلاس با خطر بسیار پایین قرار گرفت. در واقع. مقادیر پیکسل‌ها در 
نقشه نهایی به سمت صفر میل می‌کند. نتایج حاصل از کاربست دو اپراتور «یا» فازی و جمع جبری فازی 
برعکس دو اپراتور فوق‌الذکر می‌باشد. در واقع. با اعمال اپراتورهای thy‏ فازی و جمع جبری فازی» در AY‏ 
uite‏ ارزش پیکسل‌ها به سمت یک تمایل پیدا می‌کند و تقریباً کل محدوده حوضه در کلاس بسیار پرحطر 
قرار گرفت. بدین ترتیب هر چهار اپراتور مذکور برای پهنه‌بندی خطر سیلاب از کارایی پایینی برخوردار 


۱۶۰۳ نشریه جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۲ تابستان‎ Ae 


نتایج مناسب و مقبولی منجر گردید. این اپراتور ترکیبی از اپراتورهای جمع جبری فازی و ضرب فازی به حساب 
می‌آید. در عملیات GUIS‏ فازی. هنگامی که ۸ برابر ۱ باشد تلفیق مشابه جمع جبری فازی و هنگامی‌که ۸ برابر 
۰ باشد ترکیب معادل با ضرب جبری فازی می‌باشد. در پژوهش حاضر ترکیب با بهره‌گیری از گامای ۰/۷ از 
دقت بالاتری برخوردار شد. این امر با توجه به نتیجه نهایی پهنه‌بندی در سطح منطقه و مقایسه Of‏ با واقعیت‌های 
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شکل ۶- پهنه‌بندی خطر سیلاب حوضه آبریز اوجان‌چای با استفاده از منطق فازی 
Fig.4. Flood risk mapping of Ojan Chay watershed using fuzzy logic‏ 
با توجه به نقشه مذکور بالغ بر LIN‏ درصد از مساحت حوضه اوجان‌چای در کلاس با خطر بسیار 365 قرار 
گرفته است. بخش قابل‌توجهی از این پهنه‌ها منطبق بر بستر رودخانه‌های اصلی حوضه و اراضی پیرامون آن‌ها 
می‌باشد. حوضه اوجان‌چای از درجه ناهمواری بالایی برخوردار بوده و عرض دره‌ها نیز از توسعه زیادی 
برخوردار نیستند. در مقاطعی که عرض دره افزایش پیدا می‌کند دشت سیلابی در بستر آن JSS‏ گرفته که 
می‌تواند توسط سیلاب‌های مخرب تحت تأثیر قرار گیرد. کلاس با حطرپذیری زیاد نیز عمدتاً در مجاورت 
کلاس مذکور قرار گرفته و بالغ بر VV‏ درصد مساحت حوضه مطالعاتی را شامل می‌شود. این پهنه‌ها به‌طور 
عمده می‌توانند توسط سیلاب‌های با دوره بازگشت ۱۰ سال و بالاتر تحت SE‏ قرار گیرند. پهنه‌های مذکور 


عمدتاً در قسمت‌های پایین‌دست حوضه توزیع شده‌اند. افزایش عرض دره و وجود اراضی هموار در مجاورت 


سال سیزدهم اسمعیل پور ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از منطق فازی و مدل هیدرولوژیکی s‏ ۸۱ 


رودخانه‌ها از دلایل اساسی توزیع بیشتر پهنه‌های مذکور در قسمت‌های پایین‌دست می‌باشد. پهنه‌های مذکور 
در قسمت‌های میانی و بالادست حوضه در ارتباط با توسعه محلی عرض دره‌ها می‌باشد که از آسیب‌پذیری 
زیادی نسبت به سیلاب برخوردارند. بدین ترتیب مجموعه‌ای از اقدامات سازه‌ای و غیرسازه‌ای برای مهار 
سیلاب در پهنه‌های مذکور ضروری می‌باشد. در این رابطه» حدود ۳/۲ درصد از مناطق مسکونی با مساحتی 
بالغ ۱۱/۱ هکتار در کلاس با خطرپذیری بسیار زیاد و حدود 4۲/۹ درصد با مساحتی بالغ بر ۳۵/۹ هکتار در 
کلاس با خطرپذیری زیاد ably‏ شده‌اند. کلاس با خطرپذیری متوسط بالغ بر ۲۰/۳ درصد مساحت حوضه 
اوجان‌چای را شامل می‌شود. کلاس مذکور به‌طور عمده در قسمت‌های پایین‌دست حوضه توزیع شده است. 
این پهنه‌ها در اکثر موارد منطبق بر پا ی کوه‌های ارتفاعات منطقه می‌باشند و می‌توانند توسط سیلاب‌های صفحه‌ای 
(شیت‌فلادها) و در مواردی محدود توسط سیلاب‌های رودخانه‌ای تحت تأثیر قرار گیرند. خطر وقوع سیلاب 
در پهنه‌های مذکور را تا حد زیادی می‌توان از طریق اقدامات غیرسازه‌ای کاهش داد. کلاس‌های با خطرپذیری 
کم و بسیار کم به ترتیب حدود ۳۷/۳ و ۲۱ درصد مساحت حوضه را به خود اختصاص داده‌اند. کلاس‌های 
مذکور عمدتاً در قسمت‌های میانی و بالادست حوضه توزیع شده‌اند. توپوگرافی جزو مهم‌ترین عوامل مؤثر بر 
کاهش خطر وقوع سیلاب در پهنه‌های مذکور می‌باشد. ضمن اینکه در پهنه‌های مذکور سکونتگاه‌های انسانی 


محدودی وجود دارد و درنتیجه, وقوع سیلاب به مخاطره جدی تبدیل نمی‌شود. 


شکل 0- درصد مساحت هر یک از کلاس‌های be‏ سیلاب حوضه اوجان‌چای 


Fig.5. The percentage of the area of each flood risk class in the Ojan Chay watershed 
HEC-HMS شبیه‌سازی بارش - رواناب با استفاده از مدل‎ 
با شبیه‌سازی بارش - رواناب می‌توان پهنه‌های با بیشترین میزان تولید رواناب را شناسایی نمود که از نظر‎ 
اقدامات مقابله‌ای و مدیریت سیلاب مهم‌اند. شبیه‌سازی بارش- رواناب با استفاده از مدل 11۳0-11۷15 مستلزم‎ 


مدل‌سازی حوضه. هواشناسی و کنترل مشخصات دارد. در پژوهش حاضر بخش عمده‌ای از داده‌های مورد نیاز 


۱۶۰۳ نشریه جغرافیا و مخاطرات محیطی. جلد ۰۱۳ شماره ۲ تابستان‎ AY 


از طریق اکستنشن HEC-GeoHMS‏ در محیط نرم‌افزار ArcGIS‏ تهیه گردید. در این زمینه می‌توان به استخراج 
شبکه آبراهه‌ای» طولانی‌ترین آبراهه مرکز ثقل زیرحوضه‌هاء زمان تأخیر» زمان 83 GS‏ شماره منحنی ACN)‏ 
کاربری اراضی» گروه‌های هیدرولوژیک خاک» شیب زیرحوضه‌ها. شیب آپراهه‌ها؛ ارتفاع مساحت. محیط و ... 
اشاره نمود. در شکل (1) برخی از مهم‌ترین ورودی‌های داده‌های مکانی مدل ارائه شده است. لازم به ذکر است 
که توزیع مکانی برخی از ورودی‌های مهم مانند شیب و کاربری اراضی در مراحل قبلی تحقیق ارائه شده‌اند. 
در پژوهش حاضر روش‌شناسی SCS‏ به‌عنوان مبنای محاسبات بسیاری از پارامترهای مربوطه استفاده شد. 
داده‌های مورد نیاز شبیه‌سازی بارش- رواناب پس از آماده‌سازی در محیط ArcGIS‏ به نرم‌افزار HEC-HMS‏ 
ارسال گردید. ple‏ داده‌های مورد نیاز از fod‏ هیدروگراف‌ها و هایتوگراف‌ها در محیط نرم‌افزار مذکور وارد 
مدل گردید. به‌منظور شبیه‌سازی بارش- رواناب در سطح حوضه اوجان‌چای از داده‌های حداکثر بارش ۲۶ 
ساعته و همچنین هایتوگراف رگبارها استفاده شد. در جدول (۱) دوره‌های بازگشت حداکثر بارش ۲۶ ساعته 
برای حوضه آبریز اوجان‌چای SLI‏ شده است. برای محاسبه دوره‌های بازگشت از نرم‌افزار 570202 بهره گرفته 
شد. داده‌های سیلاب نیز به‌عنوان دبی مشاهداتی وارد مدل گردید. در شکل (V)‏ شماتیک حوضه آبریز اوجان- 
چای در محیط اکستنشن HEC-GeoHMS‏ نشان داده شده است. 


گروه های هیدرولوژیک 


زمان تمرکز (ساعت) 


شکل 7 مهم‌ترین ورودی‌های داده‌های مکانی مدل HEC-HMS‏ 


Fig.6. The most important spatial data inputs of HEC-HMS model 
محاسبه شده حداکثر بارش ۲۶ ساعته حوضه اوجان‌چای‎ polis -۱ جدول‎ 


Tablel- Calculated values of the maximum 24-hour rainfall of Ojan Chay watershed 


سال سیزدهم ‏ اسمعیل پون ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از منطق فازی و مدل هیدرولوژیکی  ...‏ ۸۳ 


احتمال وقوع دوره بازگشت مقدار محاسبه شده انحراف معیار 
AAY ۱۷/۰ Y ۰/0۵۰۰‏ 
/AY4 1۹/0۲ ۳ VW‏ + 
Yes ۵ Ae‏ ۰/۹۸۰ 
E ۰۹۰۰‏ ۳۳/۸۷۰ ۲/۰۹۱ 
Ya/£Y Yo ۰/۹۹۰‏ ۳/۹۰ 
Os ۰/۹/۸۰‏ 2/۷۵ 1۹۸۹ 
NN‏ ۱۰۰ 00/۸1 ۷/144 
Yre ۰/۹۹۵‏ ۹ 1/010 
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شکل ۷- شماتیک حوضه اوجان‌جای در محیط اکستنشن HEC-GeoHMS‏ همراه با نام زیرحوضه‌ها 


Fig.7. Schematic of Ojan Chay watershed in the HEC-GeoHMS extension along with the 
names of the sub-basins 


برای تخمین پارامترها در نرم‌افزار HEC-HMS‏ از فرایند بهینه‌سازی" که فرایندی تکرار شونده است. استفاده 
می‌شود. در این فرایند برای کمینه کردن میزان اختلاف هیدروگراف محاسباتی و مشاهداتی» یک تابع هدف در 
عنصری که در آن هیدروگراف مشاهداتی وجود دارد. بکار می‌رود و این تابع به‌گونه‌ای تعریف می‌شود که اگر 
مقدار هیدروگراف محاسباتی و مشاهداتی در عنصر موردنظر کاملاً برابر شوند» مقدار آن برابر صفر شود. 


همچنین» برای این تابع هدف. قیود به گونه‌ای تعریف می‌شوند که مقادیر پارامترهای تخمین زده شده در 


! - Optimization 


۱۶۰۳ شماره ۰۲ تابستان‎ Y نشریه جغرافیا و مخاطرات محیطی» جلد‎ Ag 


محدوده مجاز باشند. در این فرایند بهینه‌سازی از دو روش بهینه‌سازی یک و چندمتغیره با شروع از مقداری 
مشخص و تعیین شده توسط کاربر استفاده می‌شود و مبنای توقف آن. رسیدن مقدار تابع هدف به حد مشخص 
و يا انجام شدن تعداد تکرارهای مشخحص است )2017 .(Modaresi & Araghinejad,‏ در تحقیق حاضر از ارائه 
نتایج مدل بدون فرایند بهینه‌سازی خودداری می‌شود. در شکل (A)‏ هیدروگراف محاسباتی و مشاهداتی سیلاب 
در خروجی حوضه اوجان‌چای نمایش داده شده است. نتایج شبیه‌سازی نشان‌دهنده کارایی نسبتاً خوب مدل 
HEC-HMS‏ در شبیه‌سازی بارش- رواناب حوضه آبریز اوجان‌چای- مخصوصاً در برآورد دبی‌های پیک 
سیلاب- می‌باشد. به‌طوری که مقایسه دبی پیک هیدروگراف سیلاب مشاهداتی و محاسباتی نشان‌دهنده احتلاف 


Sink "Outlet1" Results for Run "Run 1" 


s 00 42:00 00:00 12:00 
02Dec2016 03Dec2016 
Legend (Compute Time: DATA CHANGED, RECOMPUTE) 
—— RunRun 1 ElementOutiet1 ResultObserved Flow EXPIRED ——— RunRun 1 ElementOutiet1 Result Outflow EXPIRED 
——- Run:Run 1 ElementW210 ResultOutflow EXPIRED ------ Run:Run 1 ElementR30 ResultOutflow EXPIRED 


شکل ۸- هیدروگراف خروجی مشاهداتی و محاسباتی حوضه ppl‏ اوجان‌چای 
Fig.8. Observational and computational output hydrograph of Ojan Chay watershed‏ 


سال سیزدهم اسمعیل on‏ ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از Glare‏ فازی و مدل هیدرولوژیکی ... AS‏ 


Depth (mm) 


Flow (cms) 


0.0 T T T T T T 
06:00 12:00 18:00 00:00 06:00 12:00 
02Dec2016 03Dec2016 
Legend (Compute Time: DATA CHANGED, RECOMPUTE) 
= OptOplimization 1 Element W340 ResultPrecipitation EXPIRED سس‎ OptOptimization 1 Element W340 ResultPrecipitation Loss EXPIRED 
—— OptOptimization 1 ElementW340 Result Outfiow EXPIRED ——- OptOptimization 1 ElementW340 ResultBaseflow EXPIRED 


۱ نمونه‌ای از شبیه‌سازی بارش- رواناب برای زیرحوضه شماره‎ -À شکل‎ 
Fig.9. An example of rainfall-runoff simulation for sub-basin 1 


06:00 1200 18:00 00:00 06:00 12:00 
02Dec2016 | 03Dec2016 
Legend (Compute Time: DATA CHANGED, RECOMPUTE) 
سس‎ Opt Optimization 1 Element W190 ResultPrecipitation EXPIRED m Opt Optimization 1 Element W190 ResultPrecipitation Loss EXPIRED 
سس‎ OptOptimization 1 Element W190 ResultOutflow EXPIRED — —- OptOptimization 1 Element W190 ResultBaseflow EXPIRED 


شکل ۱۰- نمونه‌ای از شبیه‌سازی بارش - رواناب برای زیرحوضه شماره ۱۵ 


۲۱2۰10 An example of rainfall-runoff simulation for sub-basin 5 


در جدول (۲) دبی‌های پیک سیلاب‌های حوضه اوجان‌چای به ازای بارش‌های با دوره بازگشت مختلف ارائه 
گردیده است. همچنین در شکل (۱۱) دبی پیک سیلاب‌های زیرحوضه‌ها به ازای بارش با دوره ناز گشت E‏ 
ساله ارائه شده است. نتایج حاصله را می‌توان به‌صورت زیر خلاصه نمود: 

- تغبیرپذیری زیادی در دبی پیک سیلاب زیرحوضه‌ها وجود دارد. در این رابطه عوامل متعددی از ded‏ 
ویژگی‌های مورفولوژیکی. پوشش زمین» شماره منحنی و ... در تغییرات شدید دبی اوج سیلاب‌های زیرحوضه- 
ها موثرند. بیشترین ميزان دبی‌های پیک سیلاب‌ها مربوط به زیرحوضه‌های شماره 9 ۲ و ۷ می‌باشد. در cpl‏ 
رابطه» دبی پیک سیلاب برای بارش با دوره بازگشت ۱۰ ساله در زیرحوضه‌های مذکور به ترتیب بالغ بر ۰/۷ 


۵ و ۷/۱ مترمکعب در انیه می‌باشد. در زیرحوضه‌های مذکور علاوه بر مساحت obj‏ عوامل متعدد دیگری 


۱۶۰۳ شماره ۲ تابستان‎ OW نشریه جغرافیا و مخاطرات محیطی. جلد‎ A 


از قبیل شماره منحنی بالا (به دلیل وجود خاک‌های با نفوذپذیری اندک و پوشش زمین) باعث افزایش رواناب- 
های ناشی از بارش می‌شود. 

- با اینکه مساحت در میزان دبی‌های پیک سیلاب‌های زیرحوضه‌ها موثر است؛ بااین‌حال همبستگی قوی بین 
مساحت و مقدار دبی پیک مشاهده نمی‌شود. این امر بیانگر این است که نقش ple‏ فاکتورهای حوضه‌ای در 
میزان تولید و نحوه انتقال رواناب از اهمیت بالایی برخوردار است. در این زمینه. شیب. کاربری اراضی» پوشش 
گیاهی و شکل زیرحوضه‌ها foe‏ تأثیر زیاد بر پارامترهایی نظیر نفوذپذیری و OW‏ بارش, زمان تخیر و 
زمان تمرکز از اهمیت بالایی برخوردارند. 

- در پهنه‌هایی از قسمت‌های میانی و بالادست حوضه به دلیل عوامل متعددی از قبیل وجود سازندهای 
زمین‌شناسی با نفوذپذیری اندک» پوشش ضعیف زمین, درجه ناهمواری VG‏ و شیب obj‏ پتانسیل تولید رواناب 
MI Sk‏ هش ان Ga Nn ESSEC OCR AR‏ 
برخوردارند. بدین ترتیب» در زیرحوضه‌های مذ کور اقدامات سازه‌ای و غیرسازه‌ای در رابطه با مهار سیلاب در 
سرچشمه‌ها لازم به نظر می‌رسد. 

- بارش‌های با دوره بازگشت بیش از ۱۰ سال می توانند منجر به سیلاب‌هایی شوند که همراه با مخاطرات 
می‌باشند. در این رابطه دبی پیک حوضه آبریز اوجان‌چای برای بارش با دوره بازگشت ۱۰ سال حدود 0٤/۷‏ 
مترمکعب در ثانیه» برای بارش با دوره بازگشت ۲۵ Sle‏ حدود 12/۲ مترمکعب در ثانیه برای بارش با دوره 
بازگشت ۵٩۰‏ ساله حدود ٩۱/٤‏ مترمکعب در انیه» برای بارش با دوره بازگشت ۱۰۰ ساله حدود ۱۲۸/۶ 
مترمکعب در ثانیه و برای بارش با دوره بازگشت ۲۰۰ ساله حدود ۱۵۲/۷ مترمکعب در ثانیه است. سیلاب‌های 
مذکور می‌توانند شهر بستان‌آباد را نیز تحت تأثیر قرار دهند. 

جدول ۲- دبی پیک سیل‌های حوضه اوجان‌چای (مترمکعب در ثانیه) برای بارش‌های با دوره‌های باز گشت مختلف 


Table2-Peak discharge of floods in Ojan Chay basin (cubic meters per second) for rains 
with different return periods 


زیر Y M T‏ سال Y‏ سال ۵ سال ۰ سال ۵ سال ۰ سال Slates‏ ۲۰۰ سال 
شماره ۱ Yo/4* ۳۱/۷۰ Yo/YY ۱/۳۹۹ ۸/40 Y/* VEM "INO‏ 
شماره ۲ ۰/۳۸ ۷ ۱/۳۷ MA‏ ۳/۹ ۳۳ ۷/5۸ 4/41 
شماره Yo/Y ۱۳/۷۳ ۸/۹۱ vié 0/4۲ ۱/۷۳ ivy iY Y‏ 
شماره ۶ 00/* ۸۱ 04/\ 1:۷ 1۷۱ "Vo‏ ۹/۰7 ۱/۸ 
شماره ۵ 1۷ EZAN ۱۳/۳۸ MOY 0/4۲ os ۱/۷۳ UN‏ 
شماره MM Y/Y AW 0/4۷ oNÉ ۱/۷۲ ۷۵ 1۵ ٦‏ 


EE ۱۳/۸۷ 4/44 VAY 0/4 YAY 1/1۲ wy ۷ شماره‎ 


سال سیزدهم اسمعیل پور ارزیابی خطر سیلاب در حوضه آبریز اوجان‌چای با استفاده از منطق فازی و مدل هیدرولوژیکی .. ۸۷ 


NEN 7/۱۱ 0/0 ¥/00 EY ۱/۳۹ 1/07 YY ^ شماره‎ 
۷/۳۹ vit £t ۳/۸ AV ۰/۸۸ AES NY ٩ شماره‎ 
۳/۰ WV VA ۱/۳۱ ۰۷۵ ۰/۱ VAR! ۰/۲ ۱۰ شماره‎ 
۱9/۸ ۱۳/۳۹۶ ۹/۷ "VAA ۷/0۱ vEt ۱/۳۹ ۹۰ ۱۱ شماره‎ 
1/7 ۱/۸ ۱۳/۶ MAY AN" Y/A Vtt ۰۹۲ ۱۲ شماره‎ 
۱/۳۰ ۸/۷۵ Vo £/0V ۹ ۱/۳۹۹ ۰/۸۰ */of ۱۳ شماره‎ 
1۱۰ YAY ZE ۱/۳۷ YA JY VAR! ۰/۸ ۱٤ شماره‎ 
A/0\ ۷/۲ 0/۰۱ EN vv) ۱/۳۱ J£ ۰۷ ۱۵ شماره‎ 
VYA ۳۳ Y/NN YNA ۱/6۰ ۹۵ ۰/۸ ۳۳ VV شماره‎ 
۸ 1۹ ۳/۹۱ ۱۳۹۷ AE ۳ UM ۰/۹ ۱۷ شماره‎ 
YoY/ YYA/YA ۹/۷ E/E 0/84 ۱۹/۹۲ Q/YA 5۷۳ خروجی‎ 
۲ 330E mM m wave m m 
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46°330°E 46°360°E 46°390°E 46°420"E 46°450°E 46°480°E 
دبی پیک‎ 
1۷ ۱:۳ (مترمکعب در انیه)‎ 


شکل ۱۱- دبی پیک زیرحوضه‌های اوجان‌جای برای بارش با دوره بازگشت ۱۰ ساله (حدود ۳۶ میلی‌متر) 


Fig.11. Peak discharge of Ojan Chay sub-basins for rainfall with a return period of 10 
years (about 34 mm) 


نتیجه گیری 
در پژوهش حاضر به‌منظور تهیه نقشه سیلاب حوضه اوجان‌چای تعداد VV‏ متغیر مکانی موردتوجه قرار گرفته 
و با کاربست منطق فازی در بستر سیستم اطلاعات جغرافیایی ترکیب شدند. نتایج نشان داد که بالغ بر ۳/۲ 


درصد از مناطق مسکونی با مساحتی در حدود ۱ هکتار در کلاس با خطرپذیری بسیار 5 2b‏ و حدود £Y/A‏ 


درصد L‏ مساحتی بالغ بر ۳0/۹ هکتار در کلاس L‏ خطرپذیری زیاد als‏ شده‌اند. حوضه اوجان‌جای دارای 


۸۸ نشریه جغرافیا و مخاطرات محیطی» جلد Y‏ شماره Y‏ تابستان ۱۶۰۳ 


توپوگرافی و شیب تندی بوده که این امر باعث محدودیت شدید فضا جهت توسعه سکونتگاه‌ها شده است. 
بدین ترتیب» بخش عمده‌ای از توسعه فضایی- کالبدی سکونتگاه‌های موجود در منطقه در بستر دره‌های اصلی 
حوضه صورت گرفته است. در سال‌های اخیر ساخت ویلاها و خانه‌باغ‌ها به پدیده‌ای qul)‏ در سطح منطقه 
تبدیل شده است. هر چند سیلاب‌های با دوره بازگشت پایین خطر جدی برای سکونتگاه‌ها و فعالیت‌های 
انسانی موجود در بستر دره‌ها ایجاد نمی‌کند؛ بااین‌حال» نتایج شبیه‌سازی بارش- رواناب با استفاده از مدل 
HEC-HMS‏ نشان داد که سیلاب‌های با دوزه باز گشت pty‏ از :۱ سال و مخضوصاً بالاتر از ۲۵ سال از uan‏ 
قابل توجهی برخوردار می‌باشند که می‌توانند بخش‌های قابل توجهی از بستر دره‌های کم‌عرض منطقه را غرقاب 
سازند. بدین ترتیب توسعه ساخت‌وساز در بستر و حریم رودخانه‌های حوضه باعث افزایش جدی خطر سیلاب 
شده است. همچنین ارزیابی متغیرهای موثر بر سیلاب و نتایج شبیه‌سازی بارش- رواناب با استفاده از مدل 
هیدرولوژیکی HEC-HMS‏ نشان داد که حوضه آبریز اوجان‌چای دارای پتانسیل تولید رواناب بالایی می‌باشد. 
نتایج تحقیق Sly‏ کارایی بالای مدل HEC-HMS‏ به‌منظور شبیه‌سازی بارش - رواناب در سطح حوضه اوجان- 
cot ua cle‏ مخصوصاً دقت این مدل در برآورد دبی‌های یک حوضه ماسب می‌باشد, به‌طوری که cis‏ 
بین دبی پیک هیدروگراف سیلاب مشاهداتی و محاسباتی کمتر از ۷ درصد می‌باشد. همچنین نتایج شبیه‌سازی 
نشان داد که هیدروگراف سیلاب زیرحوضه‌ها و حروجی حوضه اوجان‌چای تا حد زیادی تحت تأثیر ویژگی - 
های ژئومورفومتری و پوشش زمین می‌باشد. پیک‌های VE‏ در زیرحوضه‌های با کمترین میزان پوشش حفاظتی» 
شیب obj‏ نفوذپذیری کم. فراوانی برون‌زدهای سنگی و سطوح غیرقابل نفوذ مشاهده می‌شود. این زیر حوضه‌ها 
مخصوصاً در قسمت‌های نالا دست حوضه واقع شده‌اند. )2018( Yousefi, Golshan, & Pirnia,‏ کارایی مدل 
هیدرولوژیکی HEC-HMS‏ را جهت شبیه‌سازی سیل در دو حوضه کسیلیان L)‏ اقلیم مرطوب) و حوضه کارده 
(با اقلیم خشک) مورد بررسی قرار دادند. نتایج کارایی مدل نشان داد که  HEC-HMS‏ هر دو iam‏ 
آبخیز از دقت مناسبی برخوردار می‌باشد. Sy‏ کارایی مدل مذکور جهت مدیریت وقایع در حوضة peel‏ 
کسیلیان نسبت به حوضه کارده سیل نسبتاً بالاتر است. )2018( Rahimzadeh, & Habibi,‏ با استفاده از مدل 
HEC-HMS‏ حجم و دبی اوج سیلاب حوضه روانسر کرمانشاه را شبیه‌سازی کردند. ایشان نتیجه‌گیری کردند 
که درصورتی که پارامترهای متغیر مدل مانند شماره منحنی» تلفات اولیه و زمان تأخیر به‌درستی محاسبه شده 
باشند. مدل از توان بالایی جهت برآورد دبی اوج و حجم سیلاب برخوردار می‌باشد. همچنین نتایج مطالعه 
Mazidi, & Kooshki, (2015)‏ در مورد شبیه‌سازی بارش-رواناب و تخمین سیل در حوضه آبریز خرم‌آباد با 
استفاده از مدل HEC-HMS‏ نشان داد که بین هیدروگراف‌های مشاهداتی و شبیه‌سازی‌شده سازگاری وجود 
داز همست aba olas aee Doe euge‏ تاک y cb Lue‏ ریب وار اتی کر abis‏ ای مال 


را در برآورد دبی پیک در حوضه موردمطالعه نشان می‌دهد. 
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